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Abstract 

We present a partial-wave analysis of the polarization data for the reaction pp PP'^°, 
based solely on the recent measurements at lUCF for this channel. The ht leads to a 
per degree of freedom of 1.7. Methods for an improved analysis are discussed. We compare 
the extracted values to those from a meson exchange model. 


1 Introduction 

Understanding pion prodnction in nucleon-nucleon [NN) collisions near threshold is of high 
theoretical interest for various reasons. As the hrst strong inelasticity for the NN system, 
its phenomenology is closely linked to that of elastic NN scattering (for recent reviews on 
the subject of near-threshold pion production see Refs. inEi). In addition, as the pion is a 
Goldstone-boson of the chiral symmetry of strong interactions, its dynamics is strongly con¬ 
strained by this symmetry (see Ref. [H] and references therein). Recently a scheme was discussed 
that is said to lead to a convergent effective held theory even for large momentum transfer re¬ 
actions such as NN —>■ NNn H U- Conhrmation of this claim is the precondition for a 
successful analysis of the isospin violating pion production reactions measured recently, namely 
the forward-backward asymmetry in pn —>■ jB] and the total cross section measurement for 

dd —> a7r° [ 7 ]. 

A complete set of polarization observables for the reaction pp ppn^ was measured for 
the hrst time in 2001 jS]. Of the two existing advanced models of pion production in NN 
collisions pcni that include higher partial waves and therefore allow predictions for polarization 
observables, only the model of the Jiilich group um has been thoroughly confronted with 
those data. Thereby it turned out that this model failed to provide an overall satisfactory 
reproduction of these polarization observables 01121 . On the other hand, the (less complete) 
data for pp —> pmi^ as well as those for pp —> dir^ ^1] were described very well by the same 

model. So far the reason(s) for the short-coming of this phenomenological model to describe 
the neutral pion production - while being rather successful for the charged pions - is not yet 
understood^. 

^Note, however, that effective field theory studies revealed many conceptual problems in this approach, as 
discussed in |2]; it is up to now unclear how much impact those have on the description of the observables. 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 
11 
12 


Type Our notation 
T'^ 


Ss 

Ps 

Pp 


7^0 

-‘ 0 ( 00 ) 0;11 

rpO 

-' 0 ( 11 ) 0;00 

T2 

-‘ 0 ( 11 ) 2;20 

rpO 

T2 

T2 

^ 1 ( 11 ) 2;11 
T2 

T2 

^ 1 ( 11 ) 2;31 
rpl 

1 ( 11 ) 0;11 

rpl 

^ 1(11)2;11 
rn3 

1 ( 11 ) 2;31 


Notation of 

2si+lLj - 


Meyer et al. |H] 

2sf + l T . 1 
^pj ; ‘■g 


^Po 

^So 

^Do 


"Pr 


^P 
^P 

'Pi 

'Pi 

'Pi 

'P 3 


'Po,s 
'P2, S 
'Pl,p 
'Pl,p 
'P2,P 
'Pl,P 
'P2,P 
'Po,P 
'Pl,P 
'P2,P 
'P2,P 


13 

Sd 

7^2 

-^2(00)0;ll 

'A 


14 


7^2 

-‘2(00)0;31 

'R 2 

^^So,d 

15 

Ds 

T2 

-‘0(20)2;11 

'P 2 

^D2,S 

16 


T2 

-‘0(20)2;31 

^F2 

^ '^2,S 


Table 1: Partial-wave amplitudes that could contribute to pp —> pp'K° near threshold. Only 
contributions arising from the hrst 12 amplitudes were considered in the present analysis. 


The presence of spin leads to contributions of many partial-wave amplitudes, even close to 
threshold, which is the regime of interest here. It is thus difficult to draw any more concrete 
conclusion from a comparison of the model results directly with the data. It is well known 
that a partial-wave analysis is an important intermediate step towards an understanding of 
hadronic reactions: being in principle equivalent to the full data set, the partial-wave ampli¬ 
tudes can be much more easily interpreted in terms of their physics content. As a consequence, 
a comparison of the theoretical results with the partial-wave amplitudes is expected to reveal 
the strengths/weaknesses of the theory much more clearly than a direct comparison with the 
data. 

In this paper we present a first step towards a full partial-wave decomposition of the reaction 
pp —>■ ppiT^. In our work we use as input only data from the recent lUCF measurement |H]. 
However, as will be stressed below, a combined analysis of both the production data and the 
data on elastic scattering is mandatory for the future. In [S], the various angular-dependent 
structures of the polarization observables were fitted under particular assumptions on the partial 
wave content of the data as well as on the energy dependence of some of the amplitudes. These 
assumptions were necessitated by the limited statistical accuracy of the data. As we use the 
extracted coefficients of Ref. jH] as input for our fitting procedure, we also have to make the 
same assumptions in our analysis. 

This paper is organized as follows: in the next section we will describe the theoretical 
formalism that allows to relate the observables to the partial-wave amplitudes. In section 3 the 
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method of extraction as well as that for determining the uncertainties are explained. Then, in 
section 4 we discuss the results and compare them to those of a microscopic model [H]. The 
paper closes with a short summary and a discussion of further steps. 


2 Theoretical formalism 


The T-matrix for pp —> pp7r° may be expressed in the form 


1 

T= E (S\sj,s,)-T\s,,st)), (1) 

Sf,Si=0 A=|si-sj-| 


where Sj, Sf denote the initial and hnal channel-spins respectively. We use the same notations 
as in cni, where the irreducible channel-spin transition operators S^^{sf,Si) of rank A are 
dehned. If in c.m., pi,p denote the relative momenta of the two protons in the initial and hnal 
states and q the momentum of the pion, the irreducible tensor reaction-amplitudes T^^{sf,Si) 
in m can be expressed in the form HSI 


Lp,L,lgj,J,Lf * 

X TiliLps,)r,Ls.mpip) ® Yi^iq))^f 0 



Lp j 1 

J Lf] 


( 2 ) 


to separate the energy and angular dependence of the amplitudes. In Eq. we use the 
short-hand notation [j] = y 2 j^FT and (Ti 0 T 2 )^ indicates the coupling of the two irreducible 
tensors Ti and T 2 to total angular momentum L with projection m. The partial-wave amplitudes 
'^iq(LpSf)j Lsi functions of both the c.m. energy Ecm and e, the relative kinetic energy of the 
nucleon pair in the hnal state (in contrast to a two body reaction, where the partial-wave 
amplitudes are characterized by a single energy variable). 

If P, Q denote respectively the beam and target polarizations, the diherential cross section 
in a double-polarized experiment may be written as na 


da 

dhlpdrigde 


.. 1 ki+k2 

J E E ■ B\kpk^)), 

k\^k2=0 k=\k\ — k2\ 


(3) 


in terms of the irreducible tensors 


E(2v+1) E 

Si , s '.=0 A,A' 


Sf=0 


X 


r , A 1 


{s', S, fcl 

1 1 1 1 

\A A' Sf] 

1 2 2 ( 

\ 

[ fcl ^2 k) 


+A' 


(4) 
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Figure 1: The observables Aj:{0g) and Azz{0q) as a function of the pion angle at several bom¬ 
barding energies. The data and the nomenclature for the observables are taken from Meyer et 
al. jH]. The solid lines represent our results. 


which are bilinear in the irreducible tensor amplitudes whose complex conjugates 

T4(s/,iii)- define s.)'. If 

<ro«) = jBS(O.O) (5) 

denotes the unpolarized differential cross section with ^ collectively standing for {p,q,e}, the 
B^{ki,k 2 ) are related to the independent (Cartesian) spin-observables dehned in jH], 

through 



(6a) 

^..(0 = i 1) - 5i(i, 1))]; 

(6b) 

= + 72^0(1,1)]; 

(6c) 

ao(0 ^..(0 = 1) - 1)]; 

(6d) 

ao(0 AA(0 = i3fJ(i?2'(i,i)); 

(6e) 

MOAzo{0 = iB]{l,0)- 

(6f) 


(6g) 


3 Extraction of partial-wave amplitudes 

A priori, a set of 16 partial-wave amplitudes can be expected to contribute to the reaction. 
We list them in Table ^ explicitly both in our notation and the notation of Meyer et al. |H| . 
But we consider only contributions from the hrst 12 amplitudes since hnal states with orbital 
angular momentum greater than 1 were ignored in the analysis of Meyer et al. [Hj. Thus there 
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Figure 2: The observables Aj:{9p) and Azz{9p) as a function of the proton angle at several 
bombarding energies. 


are 24 real unknowns (12 complex amplitudes) to be determined. However, as overall phases 
are unobservable and s/ = 0 and s/ = 1 iViV-hnal-states do not mix with each other in any 
of the spin-observables measured in jS] (hnal state polarizations were not measured), we have 
the freedom to choose the hrst two amplitudes to be real. This leaves 22 real numbers to be 
determined. Equations (lla)-(llh) of jS] represent the general angular dependence of cto(0 and 
cro(OAi(0 in terms of the real coefhcients E, Wj, /, K, and The quantities 

E, Ek, Gfc and denote the weighted sums of bilinears in the partial-wave amplitudes cor¬ 
responding respectively to (Ss)^, (Ps)^, (PsPp) and (Pp)^ interference terms, while J, JP and 
K, FP represent respectively the contribution of (SsSd) and (SsDs) interference terms, which 
were ignored in the analysis of jS] and therefore also here. Using Eqs. (ED and (ED, we obtain 
explicit expressions for all the observables in (j(ia|l -(|bgD including the unpolarized differential 
cross section, dehned in in terms of the hrst 12 partial-wave amplitudes listed in Table ^ 
These expressions when compared with Eqs. (lla)-(llh) of [H], allow us to obtain explicitly 
the partial-wave decomposition of the coefhcients E, Ek, and . 

Note that all values given in Table IV of jS] for the various coefhcients are integrated with 
respect to the outgoing two nucleon energy, e. Thus, they are expressed as weighted sums of 
numerous bilinear in the partial-wave amplitudes, as (we use non-relativistic kinematics) 

^^max 

-BKK'(-^cm) = / T,^{Ecra,e)Tl,{Ecra,e)q{Ecra,e)p{e)de , (7) 

Jo 

where T^, k, = 2,... ,12, denotes the partial-wave amplitude listed in Table [T] (for example, 
T5{Ecm,e) = Fqipi.ii), p{e) = y/M^e and q{Ecm,e) = ^/2|J,{Ecm - 2 ,Mm - - e), with the 

reduced mass of the outgoing three body system p = 2mT,M]^/{niT, -|- 2Mn), where and 
rriTr are the nucleon and pion masses, respectively. Thus, to proceed further we need to make 
an assumption regarding the e-dependence of the T^. In the present hrst analysis of the lUCF 
data we use the most naive ansatz possible: we assume that the entire energy dependence of 
the amplitudes stems solely from the centrifugal barrier. This gives 

T,(F,^, e) oc q{E,^, e)hWp(g)UU) ^ (g) 
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Figure 3: The observables A^Q^Og), Axl{6q) and A‘^{9g) as a function of the pion angle at 
several bombarding energies. 


which should hold as long as the outgoing momenta are small compared to the inverse of the 
production radius HHI and the effects of the hnal state interaction (FSI) are negligible. (This 
is obviously wrong for the NN F-waves; they are discussed separately below.) Note, the same 
assumption was also used in the htting procedure of Ref. [S] in order to determine some of the 
coefficients E, and from the data, for the statistical accuracy of the data did not 

allow for a separate £t of these coefficients at each energy. From the ansatz of Eq. (jHl), one 
easily derives 

= z^zl, . (g) 

Thus, we find the energy dependence of the k, k' = 2,..., 12, to be of the form rj^, with 
X equal to 6,7 and 8 for the PsPs, PsPp and PpPp interference terms, respectively. Here 
h = T’max/’^TT with g^ax being the largest possible value of pion momentum for a given incident 
energy. By assumption, k = 2,..., 12 in (jH)) are energy-independent complex quantities to 
be determined from the data. 

Since the transition amplitude with the Ss hnal-state does not interfere with any of the other 
partial waves and since its FSI does not show a power law behavior na, we parameterize it as 

Bii(Ecm) = / |Ti(Ecm,e)|^g(Ecm,e)p(e)de = (10) 

Jo 

and extract it at each of the four bombarding energies individually; is directly pro¬ 

portional to the bilinear coefficient E in Table IV of Ref. [S]. 

The values of the coefficients 77)*°, F 2 , Hf, Gf’, Gy, Hf', i7|°, Hf, Hf, and 

w were determined by Meyer et ah jS] by assuming their energy dependence to be of the form 
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Figure 4: The observables Ay^{6p), Atl{6p) and A^/^{6p) as a function of the proton angle at 
several bombarding energies. 


given in Eq. ©• Therefore, in order to be consistent with jH], we do the following. We plot 
each of these 13 coefficients as a function of their appropriate r^-dependence and extract the 
slope and the corresponding error, using the values and errors in Table IV of [H]. For example, 
(<7'tot(7i)F2(?7))/(87r^) is plotted against rf. The values and errors so obtained were then used in 
our fitting procedure. This gives us a set of 13 equations for the amplitudes k = 2, ..., 12. 
The remaining 26 coefficients, viz., Fi, and the 23 coefficients from to H^, in 

Table IV of jH] lead to 26 x 4 = 104 equations for the amplitudes z^, n = 2,... , 12, as they were 
extracted by jS] without any assumption about their energy dependence. Both the values and 
errors for these 104 coefficients, taken from Table IV of jS], were multiplied by (atotiv)/ 
for consistency. In total, we have 117 equations with 21 real unknowns (since Z 2 is assumed to 
be real). In addition, the coefficient E depends only on zi. Since this amplitude is assumed 
to be real and we know the explicit dependence of E on zi, it is directly determined (up to a 
sign). The uncertainty in zi, is determined by that in E. 

This nonlinear overdetermined system of 117 equations can only have approximate solutions, 
which were obtained by y^-minimization using the software Mathematica. The resulting per 
degree of freedom was 1.7. This value was the best that we could obtain after using all the 
four methods of minimization available with Mathematica, viz., differential evolution, Nelder- 
Mead, random search and simulated annealing. We performed various checks (different starting 
vectors; setting individual amplitudes to zero) to further support that the minimum is indeed 
a total minimum. In Figs. fTlini to illustrate the quality of the fit, we compare our results to 
some of the observables measured in [S]. 

The uncertainties in the z^, = 2, , 12, were determined as follows. Let a denote the 

vector whose 21 components are the real and imaginary parts of the amplitudes Zi. Then the 
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Figure 5: Additional observables at a bombarding energy of 375 MeV as a function of the pion 
angle. 


uncertainty (standard deviation) in any of the components, say Oj, denoted by Ua., is obtained 
through 



( 11 ) 


Here fi{a) stands for the explicit functional form of the bilinear coefficients listed in Table IV 
of [H], in terms of the cxj are the corresponding errors in these bilinear coefficients, taken 
from the same Table and eji is the (j, /)-th element of the error (covariance) matrix e dehned 
as the inverse of the curvature matrix ol whose elements are given by 


_ 1 

2 dttjdai 


( 12 ) 


Here ao is the value of a for which the value of is at its minimum. 


4 Results and discussion 

In Figs. I7IIH1 the values for the z^, k = 1,..., 12, as determined in the £t are plotted. The 
uncertainties quoted in these hgures for \zk\, k = 1,..., 12 and tan(Arg(2;K)), k = 3,..., 12, 
were determined by propagating the errors obtained for the real and imaginary parts of the z^,, 
K = 1,..., 12, in the standard way [21]. It is striking that the amplitude Z 2 -, corresponding to 
the transition ^Sq — Pqs, is the largest. 
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Figure 6: Observables that depend on plotted as a function of the angle 9q at the bom¬ 
barding energies of 375 and 400 MeV, at which the measurement of Ref. [H] has the best 
statistics. 


We now turn to a comparison of the extracted with the predictions of the microscopic 
model of the Jiilich group. For a detailed description of the model we refer the reader to Refs. 0 
in]. Here we only want to summarize its salient features. In the Julich model all standard pion- 
production mechanisms (direct production (Fig. pion rescattering (Fig. IHb), contributions 
from pair diagrams (Fig. Et)) are considered. In addition, production mechanisms involving 
the excitation of the A(1232) resonance (cf. Fig. Eli-g) are taken into account explicitly. All 
NN partial waves up to orbital angular momenta Lp = 2, and all states with relative orbital 
angular momentum Iq < 2 between the NN system and the pion are considered in the hnal 
state. Furthermore all ttN partial waves up to orbital angular momenta = 1 are included 
in calculating the rescattering diagrams in Fig. |ni3,e and g. Thus, this model includes not only 
s-wave pion rescattering but also contributions from p-wave rescattering. 

The reaction NN —> NNtt is treated in a distorted wave born approximation, in the stan¬ 
dard fashion. The actual calculations are carried out in momentum space. For the distortions 
in the initial and hnal NN states, a coupled channel {NN, NA, AA) model is employed [221 
that treats the nucleon and the A degrees of freedom on equal footing. Thus, the iViV iVA 
transition amplitudes and the NN T-matrices that enter in the evaluation of the pion pro- 
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Figure 7; Left panel: Moduli (in (/ib)^/^) of Zi for bombarding energies of 325, 350, 375, and 
400 MeV. Right panel; Moduli (in (pb)^/^) of the k = 2,..., 12. 
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Figure 8: Tangent of the arguments of the z^,, k = 3,..., 12. 


duction diagrams in Fig. IHlare consistent solutions of the same (coupled-channel) Lippmann- 
Schwinger-like equation. 

By taking the partial-wave amplitudes predicted by the model, it is straightfor¬ 

ward to extract the moduli of the z^ from the model through their definition in Eqs. © and 

0 : 

Since the phases of the bilinears calculated from the model are not in all cases consistent 
with the factorization used in Eq. Q, they can not be compared easily with those extracted 
from the data. 

The \zt^\ predicted by the model are compared with the results of the partial-wave analysis 
in Fig. [ini In the upper part of the graph we compare the moduli of the of the model with 
those obtained by our partial-wave analysis, while in the lower part, the deviation of the model 
predictions from the analysis are presented. Note that the model results are normalized in 
such a way that \zi\ for a bombarding energy of 375 MeV (i.e. the one corresponding to the 
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Figure 9: Pion production mechanisms taken into account in the model of Ref. jH]: (a) direct 
production; (b) pion rescattering; (c) contributions from pair diagrams; (d) to (g) production 
involving the excitation of the A(1232) resonance, depicted by the double line. In the diagrams 
the pion (nucleon) is shown as a dashed (solid) line. 


^Pq — SqS transition) coincides with the corresponding extracted value. This is done in order 
to facilitate the comparison between the various other \zt^\. 

It is evident from Fig.^Jthat the microscopic model of Refs. jSlE] yields a rather impressive 
overall description of the various partial-wave amplitudes, cf. the hlled squares and circles. This 
is particularly remarkable because one has to keep in mind that practically all parameters of 
the model were hxed by other reactions (elastic-AA-scattering and tiN scattering)^. In fact 
the majority of the partial-wave amplitudes is reproduced even quantitatively (if one takes into 
account the error bars of the partial-wave analysis). The only serious discrepancy occurs in 
the amplitude Zg i^Pi — Sip) and to a lesser extent also in Zi 2 {^F^, — Pov)- The reason for 
the short-coming in the model prediction for these and the connection with its dynamical 
ingredients needs to be explored in the future. 

Fig. Uni suggests also considerable deviations in Z 13 , zi^ and especially in ziq. However, these 
Zi^ correspond to partial waves with NN D-waves or pion d-waves in the hnal state whose 
contribution had been set to zero in the analysis of Meyer et al. jS] — as well as in ours—as 
already mentioned above. Thus, the predictions of the microscopic model can be seen as an 
indication that those amplitudes may not be negligible and therefore should be taken into 
account in any future analysis of the reaction pp —> ppn^. Since in the present analysis, the 
neglected contributions of the D-wave and d-wave amplitudes are presumably mimicked by 
other partial-wave amplitudes, a more complete partial-wave analysis could yield results that 
are even closer to the model prediction than for the case considered in the present paper. 

^The only free parameter was fixed to the total cross section for pp —> ppir^ at low energies im 
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Figure 10: Comparison of the extracted \zk\ (filled circles) to the corresponding quantities 
predicted by the microscopic model of Ref. P (filled squares). The opaque triangles show 
results where the A contributions of the model were omitted completely whereas the opaque 
diamonds represent results of a calculation where only the A contributions after pion emission 
were omitted. In the upper panel the results for the various \zt^\ are shown, normalized with 
respect to our extracted value for \zi\ at 375 MeV, while in the lower panel, the relative 
deviations of the model calculation from the extracted values are shown. 


While a well-founded theoretical interpretation of the obtained partial-wave amplitudes calls 
for a thorough investigation, e.g. within the framework of effective field theory, as advocated 
in Refs. Elia El, the presented analysis allows already to shed light on the role of the A (1232) 
resonance for 7r° production. The importance of the A isobar for the reaction NN —> NNtt was 
already pointed out in Ref. p. The present partial-wave analysis allows to confirm that aspect 
nicely in a quantitative and transparent way. The results of the model of Ref. pea after 
omitting contributions involving A degrees of freedom are shown by the triangles in Fig. 

The corresponding predictions clearly fall short in describing the amplitudes of the partial- 
wave analysis. In particular, even the qualitative trend in the magnitude of the amplitude is 
not reproduced. 

Further insight can be gained by taking into account only those NN ^ NA transitions that 
ocenr before the pion emission (Fig. Pl and e). The corresponding predictions for the are 
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shown by the open diamonds in Fig. For almost all partial waves this part provides the 
dominant A effect, as expected, since the energy of the incoming NN system is not too far 
away from the nominal NA threshold. The energy in the outgoing NN system, on the other 
hand, is much smaller and therefore the excitation of the A(1232) is expected to be of much 
less signihcance. However, to achieve also a quantitative agreement with the extracted |zk|, the 
A excitation in the final state (after pion emission: Fig. IHF and g)) is essential, as can be most 
clearly seen in case of Z 2 Pqs) and Z 3 {^D 2 — P 2 s). Especially, only after inclusion of 

the A in the final state the former became larger than the latter. 

In this context it is also important to note that both types of contributions, the emission of 
a real pion from a A decay—as depicted in diagrams d) and f) of Fig. |H1— and the emission of a 
virtual pion from a A that gets rescattered off the other nucleon—depicted in diagrams e) and 
g) of Fig. El— are of similar numerical significance. This should not come as a surprise, for as 
soon as the Delta-isobar is involved, the large isovector pion nucleon interaction can contribute 
to the neutral pion production ^21- This is also consistent with the fact that both these 
contributions (amongst others) contribute at next-to-leading order in the chiral expansion j3] . 

In any case, it should be clear from this discussion that the A degrees of freedom have to 
be taken into account explicitly in any model that aims at a quantitative description of the 
reaction pp ppn^ even for energies near the pion production threshold. 

5 Summary and outlook 

We have presented a partial-wave analysis of the double polarization data for the reaction pp —> 
ppiT^, measured at the lUCF jH]. Due to the limited statistical accuracy of the data, following 
the authors of Ref. [Hj, we made several assumptions about the contributing amplitudes in order 
to be able to perform the analysis. The quality of the fit is with a degree of freedom 

=1.7 not completely satisfying. This could be a consequence of the several assumptions that 
were made in the analysis. 

When compared to the results of a microscopic model j2] , the analysis made three important 
points rather explicit: (i) the A degree of freedom is important for a quantitative understanding 
of the reaction pp —> ppn^, (ii) there is especially one z^ that very strongly deviates from that 
extracted from the data, namely zg {^Pi —> ^Pop )—this will guide the search for the possible 
short-comings of the model, and (hi) the set of partial waves included in the analysis was 
possibly too limited. 

As a next major step a combined analysis of NN scattering data and data on NN —>■ NNtt 
needs to be performed. On the one hand, the pion production channels provide directly the 
inelasticities to be used for the analysis of the NN data, on the other hand, the NN elastic 
phase shifts provide the phase-motion as well as the dominant e dependence of the moduli of the 
production amplitudes. The latter connection is provided by dispersion integrals as discussed 
in detail in Refs. Ca 121- 
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